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HOW  VIOLENT  ARE  FAST  CONTROLS?1 


Thomas  I.  Seidman2 


ABSTRACT: 

Consider  a  controllable  system:  x  =  Ax  +  Bu,  with  x(0)  =  0.  Given  any  time 
T  >  0  there  is  then  a  control  operator  Cj  :  £  t-t  u(-)  giving  the  (unique)  minimum 
norm  control  such  that  x(T)  =  £.  We  show  that  ||Cr||  ~  where  7  is 

computable  from  A,B  and  K  is  the  minimal  exponent  giving  the  rank  condition. 


1This  research  was  partially  supported  under  grant  AFO SR- 82-0271.  Portions  of  this  were  done 
while  the  author  was  visiting  at  the  Systems  Research  Center  (U.  of  Md.)  with  NSF  support  under 
CDR-85-00108  and  at  the  Centre  for  Mathematical  Analysis  (Australian  National  University). 

2Department  of  Mathematics,  University  of  Maryland  Baltimore  County,  Catonsville,  MD  21228. 
(BITNET:seidman@umbc;  ARPANET:seidman@umbc3). 
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1.  Introduction 

Consider  a  linear  control  system 

(1.1)  x  =  Ax  +  Ba  ,  z(0)  =  0 

with  A,B  constant  matrices  (n  x  n  and  n  x  m,  respectively,  so  x(-)  is  JRn-valued  and 
the  control  «(•)  is  IRm-valued).  Assuming  this  is  controllable,  we  know  that  for  each 
terminal  time  T  >  0  and  each  target  £  G  IRn  there  exist  controls  u(-)  giving  x(T )  =  £ 
and,  indeed,  that  there  is  unique  such  control 

(1.2)  uopt  =  uopt{-,T,  £)  G  A2([0,T]  ->  JRm)  =:U  =  UT 

minimizing  ||u||.  Of  course  the  norm  to  be  minimized  is  that  of  U  =  L2([0,  T]  — >  JRm) 
so  ||uopt||2  gives  the  least  control  energy  needed  to  reach  the  target  £  at  time  T. 

It  is  to  be  expected  that  more  violent  control  would  be  needed  as  the  time  T 
available  becomes  shorter3.  Our  object  in  this  paper  is  to  give  a  precise  (asymptotic) 
answer  to  the  question  of  the  title.  Since  the  optimal  control  uopt  is  given  by  a  linear 
operator 

(1.3)  Cr  :  £  ^  uopt{-,T,£)  :  lRn  — >  U  =  UT  , 
the  principal  result  can  be  stated  as 

(1.4)  ||CT||  ~7r~CK+1/2)  as  T  — >  0 

where  K  is  the  minimal  exponent  giving  the  well  known  rank  condition  for  control¬ 
lability: 

(1.5)  rank  [B,AB,...,AirB]  =  n 

and  7  ^  0  is  also  computable  from  A,B. 

It  is  worth  noting  that  this  also  estimates  sensitivity  for  observation  of  the  adjoint 
problem.  If  one  can  observe  y  :=  B *z  for  z  satisfying  z  =  —A *z,  then  one  easily  sees 
that  one  recovers  the  state  through  z(T)  =  C^y(-).  This  means  that  the  uncertainty 
in  the  recovered  state  due  to  noise  or  measurement  error  e(-)  in  the  observation  is 
estimated  by  || e(*) ||  <  ||Cr||  ||e(-)||.  The  same  result  (1.4)  shows  how  sensitivity  to 
error  increases  as  the  observation  time  shrinks.  Plausibly,  we  might  anticipate  that 
the  expected  noise  energy  is  proportional  to  times  i.e., 

Exp[J  \e\2dt]  ~  <j2T,  so  Exp[\\e(-)\\]  ~  crT1/2. 

3The  uniqueness  of  uopt  and  the  linearity  of  the  map:  £  i— ►  uopt  follow  from  the  Hilbert  space 
projection  theorem  under  more  general  conditions  than  here.  From  uniqueness  it  follows  that  ||uopt|| 
is  strictly  decreasing  in  T  for  each  £  ^  0. 
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The  sensitivity  estimate  then  becomes 

(1.6)  ||expected  uncertainty  in  ^(T)||  ~  70 T~K  as  T  — >  0. 

The  formula  (2.3)  is  classical  but  it  is  interesting  to  observe  historically  that  the 
question  of  the  title  seems  to  have  been  considered  first  for  distributed  parameter 
systems4  although  it  was  posed  for  the  present  finite  dimensional  case  at  least  as  far 
back  as  1975  [3]. 

2.  Formulation 

Treatment  of  (1.1)  is  expressible  in  terms  of  the  matrix  exponential,  given  by  the 
convergent  series 

OO 

(2.1)  C'A  :=  £(**/*!) A\ 

0 

The  solution  £  of  (1)  is  then  given  by 

x{t)  —  [  e<‘t~s^AlBu(s)ds 
Jo 

so,  in  particular,  one  has  x(T )  =  Vit(-)  where 

(2.2)  V  =  VT  :  U  ->  lRn  :  u(-)  fT  e(r_s)AB  u(s)ds. 

A  standard  argument  shows  ||u(-)||u  is  minimized,  subject  to  the  condition  Vm  =  (, 
by  taking  u  €  R(V*)  whence 

uopt(-’,T,£)  —  V*w  with,  necessarily,  VV*w  =  £; 

so  one  has  Ct  :  £  uopt  given  by 

(2.3)  cT  =  v;(vrv;)-1 

where  controllability  gives,  as  we  see,  invertibility  of  VyVy  =:  W7  =  W.  We  easily 
see  that  Y*  :  Mn  — *  U  is  given  by 

(2.4)  [V*w](i)  =  [e(T-t)AB]*w  for  t  G  [0,T] 
so  the  n  x  n  matrix  W  :=VV*  is  given  by 

(2.5)  W  =  W T  =  J*  esAB [esAB]*ds. 

4 One  has  log  ||Cr||  =  0(1  /T)  (sharply)  for  the  known  infinite  dimensional  cases[4],  [1],  [2]. 
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Clearly,  W  is  self-adjoint  and  (at  least)  semidefinite  from  its  form.  We  have  the 
identity 

(2.6)  !Kpt||2  =  ||Cr£||2  =  (V*W-1£,  V*W-1  £) 

=  (W-1£)  •  (VV*W_1f)  =  (Wy^)  •  £ 

which  makes  it  clear  that  our  object  must  be  to  compute  W^1  asymptotically. 

The  key  to  our  approach  is  the  invertibility  of 

(2.7)  Q  :=  lim  T^2K+^TtWtTt, 

T— >0 

using  a  suitable  family  of  operators  T  =  IV  such  that 

(2.8)  rT  invertible  for  T^O,  r7'  =  I,o  +  0(T); 
see  (2.12),  below. 

Given  the  matrices  A,B  we  consider  the  nested  sequence  (So,  Si, . . .)  of  subspaces 
of  JRn  given  recursively  by 

(2.9)  Sk  =  Sk-i  +  R(AkB)  with  S_x  :=  {0}, 

So  =  £(B),  Si  =  £(B)  +  £(AB),... 

so  each  Sk  is  the  column  space  (range)  of  the  composite  matrix  [B,  AB, . . . ,  AfcB]. 
The  assumption  of  controllability  means  that  Sk  =  -K”  for  large  enough  K  (i.e., 
(1.5))  and  we  fix  K  as  the  minimal  exponent/index  giving  this. 

For  each  k  (0  <  k  <  K)  we  can  find  the  orthogonal  complement  of  Sk-i  in  Sk 
and  let  Efc  be  the  orthogonal  projection  on  this  subspace.  This  gives  the  important 
fact  that 

(2.10)  EiAJB  =  0  for  j  <  k  <  K 

since  j  <  k  gives  £(AJB)  C  Sk-i  C  M(Bk).  We  observe,  although  we  do  not  need 
the  fact,  that 

m  >  dim  R  (B)  =  dim  R  (Eo)  >  dim  R  (Ei)  >  . . .  >  dim  R  (E# ) ; 

we  will  need  the  fact  that  Sk-i  7 -  Sk  —  by  the  definition  of  K  so  dim  R(EK)  ±  0 
and  E k  7^  0.  The  construction  of  {E*}  gives  a  direct  sum  decomposition 

1  =  Eo  +  . . .  +  E k  ,  m  ~  ©o  ^  5 

(2.11)  Sk  —  £(Eo)  ©  ...  ©  £(Efc)  for  k  =  0,...,K. 

Thus,  introducing 

(2.12)  T  =  TT:=itk'-TK~kEk 

o 

we  see  that  (2.8)  holds  with  T0  =  K\Bk  7^ 
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3.  Principal  Computation 


Our  object  in  this  section  is  to  obtain  (2.7),  with  as  in  (2.12),  computing  Q 
and  showing  it  is  invertible. 

The  integral  expression  (2.5)  gives,  on  substituting  s  =  To, 

r-(2ic+i)rwr  _  f1[T~KTeT(7AB}[T-KTeTaAB}*d(T. 

Jo 

Using  (2.12)  and  (2.1),  we  have  (for  T  >  0) 

K  oo  j 

T~KYeTaAB  =  yk\T~kBk  y—T3A3B 
k=0  3=0  J' 

K  °°  lclrri 

=  EEt^'b. 

k= 0  3=0  J  • 


By  (2.10),  the  terms  with  j  <  k  vanish  so  no  negative  powers  of  T  actually  appear  on 
the  right;  we  then  split  the  sum  into  the  terms  with  j  —  k  and  those  with  j  >  k  +  1 
for  which  we  set  i  =  j  —  (k  +  1)  =0,1,  —  Thus, 


(3.1) 


T~KTeTaAB 


= 

k= o 

'  K  oo  K]  i+k+lrni 

+t  yy y — 7 — -e.a^b 
ItitA'  +  *+!)! 

=  PH  +  rR,(r,<r). 


Restricting  our  attention  to  T  <  1,  which  is  certainly  permissible  as  we  are  only 
interested  in  the  limit  T  — >  0,  an  easy  estimation  gives  the  uniform  bound 


K  oo 

|RiP»ll  <  EEIIA||i+‘+1l!B!l/(.'+i)! 

k=0 i= 0 

=  (i+...+iiAiniiBH(e«A»-i) 


since  (z  +  k  +  1)!  >  ( i  +  l)!  A:!.  Hence,  (3.1)  gives5 

T  kT  eTaAB  =  P(er)  +  0(T ) 


and 

(3.2)  {T~KPeTaAB)  (T_irreT<7AB)* 


=  [P (a)  +  TR^Ta)]  [P (a)  +  TR^T,*)]* 

=  P(a)P*(a)  +  TR2(T,a), 


5Note  that  our  estimation  of  Ri  precisely  legitimates  the  use  of  the  0(T)  notation. 
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with  R2(T,  a)  :=  (RiP*  +PR^  +  TR1R1)  uniformly  bounded.  Thus,  integrating, 

(3.3)  r(2if+i)rwr  =  Q  +  (5(r) 
with  0(T)  =  TfR2dcr  =:  TR3(T)  and 

(3.4)  Q  :=  (' V{o)Y*{o)do 

Jo 

K 

=  Y.  (j  +  k  +  1)“1Ej-Aj'BB*A*A:Ea. 

j,k-u 

We  must  show  that  Q  is  invertible. 

Lemma  :  B*A*A:Ea:£  =  0  Efc£  =  0. 

PROOF  :  For  any  £  6  IRn  we  have  E *£  £  Sk  :=  Sk-i  +  £(A*B)  by  definition  so  we 
may  write 

E*£  =  AkBr)  +  £'  (£'  £  Sk-i) 

for  some  r]  £  ]Rm.  Then,  assuming  B*A*A:Ej t£  =  0,  we  would  have 

IIEifU*  =  (A*Bt)  +  {*)  •  (Eif) 

=  n  ■  (B’A**E»e)  +  (B,0  •  £  =  0 

since  E*£'  =  0  for  £'  6  Sk-i-  □ 

From  (3.4)  we  see  that 

6-Q£=  /1£-[P(a)P*(a)£]dr=  f1  ||P*(o)£||2dr 
Jo  Jo 

so  Q£  =  0  only  if  P*(er)£  =  0. 

From  the  definition  of  P(-),  this  would  mean  that  each  term  crA:B*A*fcE(fe£  would 
have  to  vanish  and,  by  the  Lemma,  this  would  imply  E *£  =  0  for  each  k.  Hence, 
from  (2.11),  Q£  =  0  would  give  £  =  0.  We  have  thus  shown  that  Q£  =  0  only  for 
£  =  0;  for  an  n  x  n  matrix  Q,  this  ensures  invertibility. 

4.  Results 

We  must  draw  the  desired  conclusions  from  (3.3). 

It  is  clear  from  the  bound  on  Ri(T,  o)  and  the  obvious  fact  that  P(cr)  is  bounded 
uniformly  on  [0, 1]  that  R2(T,  cr)  is  uniformly  bounded  so  R3(T)  is  uniformly  bounded 
-  say,  ||R3(T)||  <  M3  for  0  <  T  <  1.  Restricting  attention  to  T  <  l/2M3||Q~1||  =:  r, 
we  have 


(Q  +  TR3)-1  =  Q-^l  +  TQ-'Rs)"1 
=  Q"1+TR4(T) 
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with 


IIR.II  <  iiQ-1||||(i  +  rQ-‘R3)-I-i||/r 

=  IIQ-'ll  ||Q-1R3(l  +  rQ-1R3)-1|| 

<  ||Q_1||  ||Q-1R3||/(l-r||Q-1R,||) 

<  2M3||Q_1||2  for  0  <  T  <  t, 

Now,  inverting  each  side  of  (3.3)  is  legitimate  for  T  ^  0  and  gives 

=  Q_1  +  rR4, 

(4.1)  T21 c+1Wy1  =  rr(Q-1  +TR4)rT 

=  ToQ-'To  +  TRsiT) 

with,  obviously,  ||R5(T)||  uniformly  bounded  on  0  <  T  <  t.  In  particular,  this  proves 
(independently  of  the  standard  controllability  arguments)  the  invertibility  of  Wy  — 
at  least  for  small  T  >  0  and  so  a  fortiori  for  all  T  >  0  by  the  non-negativity  of  the 
integrand  in  (2.5). 

From  (2.6)  and  the  positivity  of  W,  W-1  we  then  have 

(4.2)  ||Cr|[2  :=  max{||Cr£||2  :  ||£||  =  1} 

=  max{e-W^M|e||  =  l}  =  ||W^|| 

=  r-(2*'+i)[||r0Q-1ro||  +  o(T’)]. 

This,  of  course,  is  just  (1.4)  with,  from  (2.12), 

(4.3)  7  :=  ||roQ_1ro||1/2  =  K'.WEkQ-'-EkW112 

once  one  show  E^-Q_1E/<-  ^  0  so  7  /  0.  Note  that  the  positivity  of  Q,  hence  of  Q-1, 
gives 

||EjfirQ-1Ejfir||  =  max{£  •  (E^Q^E^)  :  \\t\\  =  1} 

=  max{£  ■  Q-1£  :  ||f||  =  1,  £  =  Ejr£} 

=  ||Q_1U(ek)||- 

Since  j?(Ejf)  7^  {0}  by  the  minimality  of  K,  this  is  clearly  non-zero. 

At  this  point  we  work  out  in  somewhat  greater  detail  the  case  of  scalar  control 
(m  =  l) .  The  n  x  1  matrix  B  is  now  just  a  vector  and  we  set 

(4.4)  Po  —  B,  Pk  =  A*/?0  for  k  =  0, . . . ,  n  —  1. 

Note  that  controllability  gives  K  =  n  —  1  in  this  case  so,  for  scalar  control,  our  result 
(4.2)  becomes 

(4.5)  ||Cr||  ~  7 T-<n+1/2)  + 
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To  compute  7  here,  note  first  that  (/30, . . . ,  /3n_i)  is  a  basis  for  ]Rn  and  let  (er0, . . . ,  en_i) 
be  the  orthonormal  basis  obtained  from  that  by  the  Gram-Schmidt  procedure  so 
E*  :  £  (f  ‘  £k)£k  and  B*A*fcEfc£  becomes  (/?*,  ■  er*,)  (£*,  •  ek).  Then  (3.4)  becomes 


(4.6) 


K-l 

Q(=H 

j=0  L 


k= 0 


j  +  k  +  1 


This  shows  that,  re-written  in  terms  of  the  orthonormal  basis  (£0, . . .  , £n_i),  the  new 
matrix  for  Q  is  just  DHD  where  D  :=  diag  [flj  ■  Sj]  and  H  is  the  n  x  n  Hilbert 
matrix.  Then 


||EjfrQ“1EJC 


=  £K  '  Q  l£K 

—  [lower  right  corner  element  of  D-1H-1D-1] 

—  (/?„_ i  •  e„_i)2  [lower  right  corner  element  of  H_1] 


whence 

(4.7)  7  =  Cn(pn- 1  •  er„_i) 

with  Cn  :=  (n  —  1)!  [lower  right  corner  element  of  the  inverse  of  the  HR1]1/2.  The  co¬ 
efficient  Cn  grows  extremely  rapidly  with  n  but,  of  course,  is  fixed  for  any  given 
dimensionality.  Thus,  /?n_i  ■  £■„_!  provides  the  only  dependence  on  the  particu¬ 
lar  system  (1.1);  it  is  just  the  norm  of  the  component  of  An-1/?o  orthogonal  to 
span{/30,. . . ,  An-2/?0}. 

Returning  to  the  general  case,  we  now  consider  the  asymptotics  for  a  particular 
target  £  (rather  than  the  ‘worst  case’  treatment  above).  We  have,  from  (2.6), 

x(T;u(-))  =  (  =)•  ||«(-)||  >  ||Crf||  =  ((  ■  Wf'e)1/2. 


From  (4.1)  we  have 

||Cr£||  =  T-{k+1/2)K\{Zk  ■  Q-^k)1/2  +  0{T^K~1/2)) 
where  we  have  abbreviated  £k  :=  Ej<-£.  We  may  write  this,  assuming6  £k  /  0,  as 

(4.8)  IlC^eil  ~  {K\\\Q~1/2 U\\)T~{K+1/2) , 

which  gives  the  same  asymptotic  growth  rate  for  (almost  all)  targets. 

6 This  is  ‘almost  always’  true  —  it  fails  only  when  £  happens  to  lie  exactly  in  the  (proper)  subspace 
>/(Eji:)  in  which  case  one  has  slower  blowup.  Even  in  that  case  slight  perturbations  would,  almost 
inevitably,  give  some  component  in  £(E#)  so  this  analysis  would  dominate. 


9 


References 

[1]  E.N.  Guichal,  A  lower  bound  of  the  norm  of  the  control  operator  for  the  heat 
equation,  J.  Math.  Anal.  Appl.  110,  pp.  519-527  (1985). 

[2]  W.  Krabs,  G.  Leugering,  and  T.I.  Seidman,  On  boundary  controllability  of  a 
vibrating  plate,  Appl.  Math.  Opt.  13,  pp.  205-229  (1985). 

[3]  T.I.  Seidman,  Boundary  control  and  observation  for  the  heat  equation,  in  “Calcu¬ 
lus  of  Variations  and  Control  Theory”  (D.L.  Russell,  ed.),  pp.  321-351,  Academic 
Press,  N.Y.,  (1976). 

[4]  T.I.  Seidman,  Two  results  on  exact  boundary  control  of  parabolic  equations,  Appl. 
Math.  Opt.  11,  pp.  145-152  (1984). 


